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ABSTRACT 

The direct detection of binary systems in wide-field surveys is limited by the size of the stars' point- 
spread-functions (PSFs). A search for elongated objects can find closer companions, but is limited by 
the precision to which the PSF shape can be calibrated for individual stars. We have developed the 
BinaryFinder algorithm to search for close binaries by using precision measurements of PSF ellipticity 
across wide-field survey images. We show that the algorithm is capable of reliably detecting binary 
systems down to ~l/5 of the seeing limit, and can directly measure the systems' position angles, 
separations and contrast ratios. To verify the algorithm's performance we evaluated 100,000 objects 
in Palomar Transient Factory (PTE) wide-field-survey data for signs of binarity, and then used the 
Robo-AO robotic laser adaptive optics system to verify the parameters of 44 high-confidence targets. 
We show that BinaryFinder correctly predicts the presence of close companions with a <5% false- 
positive rate, measures the detected binaries' position angles within 2° and separations within 25%, 
and weakly constrains their contrast ratios. When applied to the full PTE dataset, we estimate 
that BinaryFinder will discover and characterize ^450,000 physically-associated binary systems with 
separations <2 arcseconds and magnitudes brighter than m^ = 18. New wide-field synoptic surveys 
with high sensitivity and sub-arcsecond angular resolution, such as LSST, will allow BinaryFinder to 
reliably detect millions of very faint binary systems with separations as small as 0.1 arcseconds. 
Subject headings: binaries: close; methods; data analysis; techniques: image processing; surveys; stars: 
statistics 



1. INTRODUCTION 

The development of ground-based high-angular- 
resolution ima ging tech niques such as adap tive optics 
(e.g. fWizinowich et al. 2006; Hcrriot et al. [2000; Troy 



et al.r 


20001 


Rousset 


et al.| 120001 iMoretti et al.||2009t 


Neiche 


1 et a 


.l'2010; Hart et al. 2010), Lucky Imaging 


(e.g. Law et al. 2006 


) and aperture-masking (e.g. |Mon- 


nier et al.||2004, ^Lacour et al.,|2011 


Tuthill et al. |liOOO|) 



has been driven by the need to discover or character 
ize closely-separated objects, usually after initial target 
identification by a wide-field survey. Such follow-up is 
typically limited by oversubscription and telescope effi- 
ciency to at most several hundred tar gets per survey (for 
example, in the low-mass-star regime: [Siegler et al.|2005{ 
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Close et al.]|2003| [Law et al.||2010[ [Janson et ar]|2012 ). 

New development s such as the high-eth ciency robotic 
adaptive optics system Robo-AO (Baranec et al. 20121 
are increasing the possible sample sizes to thousands of 
targets. However, there remain many orders of magni- 
tude more targets which cannot feasibly be searched one- 
by-one. 

On the other hand, synoptic wide-field imaging covers 
millions of objects in each exposure, but with a seeing- 
limited resolution of 0.5-2 arcseconds. The P alomar 
Transient F actory ([Law et al. 2009b Rau et al. 
Pan-Starrs ( Kaiser eF^Tp^™"^ 



J2009 

tSkyma pper (^Kelle r et al. 



2007), and ultimately LSHTllWcet al. 2008; LSSTEcT 



ence Collaborations et al. 2009 ) are, or will be, generating 
extremely large datasets covering billions of objects. 

In this paper we present BinaryFinder, a technique 
which is capable of efficiently searching wide-field synop- 
tic survey data for close companions. The techniq ue uses 



algorithms developed for w eak-lensing surveys ( Kaiser 



et al.||1995[ jHoeks tra et al.|[l998, 2005^ to measure the 
ellipticity of individual target star images. We extend 
those techniques to obtain the probability of each imaged 
star having a close (0.1-2.0 arcseconds) companion and 
to determine the position angle, separation and contrast 
ratio of high-confidence binary systems. 

Simple measurement of the ellipticity of stars in wide- 
field imaging data can reveal binaries of moderate sep- 
aration (e.g. Kraus fc Hillenbrand||200T ). However, im- 
perfect optics and atmospheric effects lead to changing 
point spread functions (PSFs) across a wide field (PSF 
anisotropy), preventing the easy detection of close and/or 
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high-contrast binaries. Our method uses hundreds of 
point sources in the field as cahbration objects to mea- 
sure PSF anisotropy, and multiple-epoch data to correct 
for the effects of changing seeing and to measure the sep- 
arations and contrast ratios of the detected systems. 

The paper is organized as follows. In Section [2] we de- 
scribe the PSF-ellipticity measurement and anisotropy 
correction algorithm, as well as the practical considera- 
tions of implementing the method with wide-field survey 
data. Section [3] describes our method of generating ar- 
tificial point sources and binary systems, applying sim- 
ulated PSF anisotropics, accounting for instrumental ef- 
fects, and matching the resulting ellipticity distributions 
to those seen on-sky. A verification of our method is pre- 
sented in Section |4] where we use Robo-AO to confirm 
and characterize binaries discovered by our algorithm. 
In Section [5] we explore the use of multiple-epoch data to 
measure the position angle, separation and contrast ra- 
tios of detected binaries. We conclude in Section |6] with 
a discussion of the false-positive rate of the algorithm, 
and the applications of the method. 

2. THE BINARYFINDER ALGORITHM 

The BinaryFinder algorithm can be divided into these 
steps: 

1. Measure the raw ellipticities and PSF anisotropy 
parameters of al l sta rs in a wide-field seeing-limited 



image (Sections 2.1 & 2.2) 



Divide the image into smaller subsections and bin 
the stars within each subsection into groups of sim- 
ilar flux (Section |2.3[ ). For each of these groups, de- 
rive a polynomial fit which describes how the PSFs 
vary as a function of position on the CCD, and use 
the fit to correct the raw ellipticity measurements 
to o btain anisotropy-corrected ellipticities (Section 
211). 



Repeat this process for multiple exposures of the 
same field, providing several ellipticity measure- 
ments for each object in the field; consolidate the 
individual ellipticity measurements of each object 
by finding the seeing dependence of the ellipticity 
and calculate t he el lipticity at a common reference 
seeing (Section 2.5). 



4. Evaluate each obj ect f or the presence of close com- 
panions (Sections 2.6 & 2.7 1. 



5. Use the changing ellipticity as a function of seeing 
to constrain the separations and contrast ratios of 
individual systems (detailed in Section [5]). 

In the following sections we describe each component in 
detail, with a particular application to wide-field survey 
data taken by the Palomar Transient Factory. 



2.1. Ellipticity and PSF anisotropy measurement 
Hoekstra et al. (2005) (hereafter HWU) describes a 



method to measure the ellipticity of a single object with 
corrections for PSF anisotropy, which we summarize 
here. HWU quantify the ellipticity of a PSF with el- 
lipticity parameters ei and 62, which relate the second 



moments of flux as measured from the centroid of the 
object. 



ei 



hi - I: 



22 



hi 



,62 



2/12 



'22 



hi 



'22 



where 



In =^/(x,2/)x2W^ 



^f{x,y)xyW 



(1) 



(2) 



(3) 



(4) 



Here, x and y are the x and y pixel numbers as mea- 
sured in a coordinate system whose origin is at the cen- 
troid of the object, f(x,y) is the source flux falling on 
the pixel located at the point (x,y) in that same coordi- 
nate system , and W is a weight function which accounts 
for the increase in Poisson noise with decreasing photon 
count away from the centroid. 

The ellipticity of the PSF of an unresolved binary or 
an extended object is a combined measure of the true el- 
lipticity caused by the astrophysical configuration of the 
object, and the ellipticity induced by directional smear- 
ing of starlight by the atmosphere and telescope optics. 
Thus, any raw ellipticity parameters need to be corrected 
for the induced ellipticity in order to obtain astrophys- 
ically significant measures of the true ellipticity. HWU 
do this using the smear polarizability tensor P which 
quantifies the changes in PSF ellipticity resulting from 
perturbations in the PSF anisotropy. 



Plj -^ij ^^^j 



(5) 



where 



Xn = T^-r- E /(^' + 2W^'(^W) + W"{x'-y')) 

^11 + -'22 _ 

(6) 



x,y 



^11 + ^22 1;^ 

(7) 
(8) 

eT - ^^-^ E /(^' y)(2^y)(2W^' + W'^x" y^)) 
^"+^22 ^ 

(9) 

Here, the differentiation of the weight function W is 
taken to be with respect to a? -/-i/2 _ 

The smear polarizability tensor of a source allows 
HWU to measure the PSF anisotropy parameters p\ and 
Pi which quantify the anisotropy in the image. 



P 

1.1. 



(10) 
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To correct the ellipticity of a star in an image, it is nec- 
essary to examine how the anisotropy parameters vary 
across the CCD in the vicinity of the star. HWU mea- 
sure the anisotropy parameters of point sources near the 
star of interest, and create a polynomial fit of pi and P2 
as a function of position on the CCD. They then correct 
the observed ellipticity of the star via the relationship 



PiiPi - a(pi +P2) 



2\0.5 



(11) 



where a is a constant proportional to the magnitude of 
the anisotropy at the point on the CCD where the star 
is located. This provides a final, anisotropy-corrected 
measure of the PSF ellipticity. 

HWU measure the ellipticity of the object in many sep- 
arate observations. The observations are taken at vary- 
ing seeings, thus they determine the ellipticity at a refer- 
ence FWHM by deriving a fit of ellipticity as a function 
of FWHM and evaluating the fit at the chosen reference 
radius. 

2.2. Object detection and measurement of raw PSF 
parameters 

In order to measure the necessary parameters of all 
the stars in our images (using Equations 1-10), we re- 
quire centroid coordinates and Full- Width at Half-Max 
(FWHM) measurements of the PSF of each source. We 
obtain a catalog of objects (these include stars, blended 
binary systems, foreground/background star ble nds, and 
galaxies) for each image using SExtractor (Bertin & 
Arnouts||1996). SExtractor provides the sky coordinates 



(RA and DEC) and pixel coordinates (x and y) of the 
centroid of each object, the FWHM assuming a Gaus- 
sian core, the fiux, and the SExtractor error fiag. We re- 
move all objects with non-zero SExtractor flags to avoid 
saturation, bad pixels, and the wider-separation binaries 
which SExtractor can directly detect and deblend. 

Using the catalog of objects obtained from SExtractor, 
we measure the raw ei and 62 ellipticity parameters and 
the PSF anisotropy parameters pi and p2 of every object 
in an image. The weight function we employ to account 
for changing SNR across the PSF is a Gaussian with 
dispersion equal to the FWHM calculated by SExtractor, 
centered at the SExtractor-derived centroid coordinates. 
We define an object aperture with a radius equal to twice 
the FWHM, and a sky-subtraction aperture with radii 
between 2.5x and 3.5x the FWHM. 

2.3. Image partitioning 

The anisotropy profile of images can vary significantly 
in wide-field images due to geometric distortions. In PTE 
data, we found it necessary to partition the images into 
subsections, and derive separate fits for each partition 
(See Figure [1]). It is important to make the regions as 
close to square in shape as possible, as an accurate fit 
requires a comparable span in each direction. Addition- 
ally, the subset of sources used to create the fits should 
always extend further than the boundaries of the region 
being corrected (as shown by the blue square in Figure 
[1]), to ensure that the ellipticities of all sources are cor- 
rected based on a fit obtained from sources surrounding 
them on all sides. 

2.4. PSF anisotropy fit 



Fig. 1. — Partitioning of image (part of PTF dataset). The im- 
ages are 2048 by 4096 pixels with a 1.01" pixel size. The red grid 
divides the image into 8 smaller subsections about 800 by 800 pix- 
els in size, each of which will be corrected for anisotropy with a 
separate polynomial fit. The larger (1200x1200) blue square which 
encloses one of the partitions of the red grid shows the region from 
which the fit will be derived for that particular partition. 



We next derive a polynomial fit of the anisotropy 
parameters for each group to determine how the PSF 
anisotropy varies across the CCD. We perform a Bcr clip 
on the values of ei, 62, pi and pi to remove severe 
outliers. The points that remain may not all be point 
sources (some are blends, galaxies, etc). In contrast to 
instrumentally-induced ellipticity, this will not cause any 
directional bias in the fit, as we expect astrophysically- 
extended objects to appear with equal probabilities in all 
orientations. 
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Our anisotropy fit is of the form 



fit 

PI = CiX- 



■C2X +C3y + Ciy +C5xy + CQ' 



'Cjftot+Csftot 

(12) 

where are the constant coefficients to be found which 
minimize the residuals of the measured pi parameters of 
the objects being fitted and x and y are the CCD pixel 
coordinates where t he objects are located. 

We extended the Hoekstra et al. (2005) model to in- 
clude terms which are proportional to the flux of the ob- 
jects (/fot). While applying the algorithm to the PTF im- 
ages, we noticed that the ellipticities of high-flux objects 
were not corrected accurately to zero using a fit depen- 
dent only on x and y object coordinates. Upon closer 
examination of the images, we realized that for high 
flux (yet unsaturated) objects, slightly-less-than 100% 
Charge Transfer Efficiency (CTE) leads to a smearing in 
the direction of readout. This causes an additional elon- 
gation dependent on the flux of the object, requiring flux 
terms in the anisotropy flt. 

A second-order polynomial fit appears to be sufficient 
for PTF images; first degree polynomial fits were insuf- 
ficient to accurately quantify the polarization variation, 
and we found that third degree and higher fits do not 
increase the accuracy of the corrected parameters com- 
pared to the second order fit. 

A similar fit is also derived for the p2 parameter, al- 
though the flux-dependent terms are not necessary as 
readout-induced elongation exactly parallel to the pixel 
grid (whether it is the x or y axis) does not affect the 62 
parameter, by definition (see equations 1 through 4). 

Once we have derived these fits we evaluate them at 
the coordinates of the objects to be corrected and obtain 
the corrected ellipticities according to Equation 12. For 
the purposes of evaluating the ellipticity of an object 

compared to others in the field, we use e, the \/ef + 62 
combination of the ellipticities. 

We find that the typical instrument-induced ellipticity 
has a magnitude between 0.01 and 0.1, roughly the same 
amount as our final measured ellipticities for detected 
binary systems. 

2.5. Seeing correction 

When we measure the corrected ellipticities of an ob- 
ject across multiple images, we measure many values for 
ei and 62 at a variety of FWHMs. However, ellipticity 
is not constant with seeing; close objects induce smaller 
ellipticities when the PSFs are large. To make ellipticity 
measurements comparable between fields with different 
seeing distributions, we create a fit of object ellipticities 
as a function of FWHM, and evaluate this fit at a refer- 
ence seeing value (Figure [2]) . 

2.6. Selecting suitable fields and objects 

Undersampled images with FWHMs close to the pixel 
scale of the CCD cannot accurately measure ellipticities 
because the ellipticity measurement is then very sensi- 
tive to the location of the object relative to the pixel 
grid. We found that pixel effects become negligible for 
roughly Nyquist-sampled images with at least two pixels 
across the star PSFs, and we therefore use only those 
well-sampled PTF images. We require at least 30 sep- 
arate observations of a target in order for us to obtain 




2.4 2.6 

seeing (arcsec) 



Fig. 2. — Seeing fit of ei, 62 and e of an elongated PTF object. 
Tlie points represent the various ellipticity measurements obtained 
from each unique image of the PTF field the object is in. The curve 
shows the polynomial fit of ellipticity as a function of the FWHM. 
The fits of all objects are evaluated at a common reference FWHM 
to make ellipticities comparable across fields. 

an accurate fit of ellipticity with seeing, and we also re- 
quire several hundred sufficiently bright sources within 
our image for PSF measurement (the last requirement 
is satisfied by essentially all PTF images). For R-band 
PTF images, our simulations and on-sky measurements 
showed that the faintest useful sources have m^ ~ 18. 

2.7. Results for a typical wide-field-survey image 
Figure Is] provides an example of the distribution of ei. 



62 and e (all measured at a common reference FWHM) 
of the objects appearing in a single chip of a PTF field 
(1/12 of a PTF field). The distributions have two compo- 
nents - the bulk of the objects which are closely centered 
around zero ellipticity, and the extended wings. The first 
component represents true point sources or blends of ob- 
jects too close together for our algorithm to differentiate 
from point sources. The thickness of this component of 
the distribution is due to measurement errors in the el- 
lipticities. The wings of the distribution represent either 
blends of multiple point sources, or other extended ob- 
jects such as galaxies. Figure [3] shows a sample of PTF 
objects from the same chip, ordered by increasing ellip- 
ticity. 

3. SIMULATED PERFORMANCE 

We tested the sensitivity of BinaryFinder using sim- 
ulated wide-field-survey data. For each simulated im- 
age we generated 2000 point sources, along with artifi- 
cial binary systems of varying separations and bright- 
ness ratios. We modelled the PSF of each object using 
a Moffat function. Then, we convolved all the object 
PSFs with anisotropic kernels to simulate the instrumen- 
tal PSF anisotropy. In order to make this anisotropy de- 
pendent on the location within the CCD, we placed four 
reference anisotropy kernels at the four corners of our 
CCD. We convolved each object with a kernel which was 
the weighted average of the four corner kernels, with the 
weights proportional to the inverse of the distance from 
the object centroid to that corner of the CCD. We then 
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Fig. 3. — Distributions of ei, 62 and e of objects in a PTF field. 
The ellipticity parameters have all been evaluated at a reference 
FWHM of 2.5 arcseconds 



e=0.0005 
■ 


e=0.004 

i 


e = 0.011 
1 


e=0.021 


e=0.032 

*- 

■ ■ ■ 


e = 0.041 


e=0.061 


e=0.15 


e = 0.22 










5 10 15 20 5 10 15 20 5 10 15 
X (asec) 



Fig. 4. — A sample of PTF objects ordered by increasing ellip- 
ticity. Each image is 20 by 20 arcseconds in size and has a FWHM 
of ~2.4 arcseconds. 




flux ratio 

Fig. 5. — Separation and contrast dependence of the ellipticity 
of the PSF of blended binaries (artificial images). The black curve 
represents the binary detection criterion (e=0.02). 

applied BinaryFinder, using the artificial point sources 
to derive the anisotropy fits and correct the raw elliptici- 
ties. The changes in the ellipticity measurement with the 
brightness ratio and separation of our artificial binaries is 
shown in FigurejSj As expected, the anisotropy-corrected 
ellipticity increases for binaries of greater separation and 
greater brightness ratios, and vice versa. 

We also utilized artificially generated point sources to 
determine an ellipticity boundary point between likely 
point sources and secure binary detections. We chose 
1200 PTF objects at random, and created 1200 similar 
artificial point sources using the method outlined above. 
For the artificial sources, we simulated the observation 
conditions and data reduction of the real images, includ- 
ing variations in cloud cover and background brightness 
which change the signal to noise ratio from image to im- 
age. We also introduced a per-pixel flat field noise in the 
simulated images at the 1.4% level. We plotted the ellip- 
ticities of the artificial point sources and the real objects 
(see Figure[6]) to compare the two distributions. The sim- 
ulations accurately reproduce the PTF-measured distri- 
bution (within sampling noise) in the area corresponding 
to point sources. However, the PTF objects have a much 
larger tail in ellipticity, corresponding to astrophysically- 
extended objects. We found that the frequency of sim- 
ulated point sources falls to zero for e above 0.02 in our 
simulations, and we adopt that as the boundary point 
between the point sources and the wide separation bi- 
naries. Note that the ellipticity values of the boundary 
point are pixel scale dependent and thus have to be cal- 
culated independently for new datasets. 

4. VERIFYING BINARYFINDER DISCOVERIES WITH 
ROBO-AO ADAPTIVE OPTICS IMAGING 

To verify the performance of the BinaryFinder algo- 
rithm we used the Robo-AO robotic laser adaptive op- 
tics system to obtain 0.1-arcsecond resolution images of 
a large sample of possible binaries. To generate the tar- 
get list we measured the astrophysical ellipticities of ap- 
proximately 100,000 objects in PTF data covering a w70 
square degree area (see Figure [t] for the distribution of 
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Fig. 6. — EUipticity distributions of 1200 artificial point sources 
and 1200 randomly chosen PTF objects (top - ei, middle - 62 , 
bottom - e). As can be seen in the histograms, the frequency of 
the artificial point sources goes to zero for magnitudes of el or 
e2 greater than KiO.015 or for e greater than KiO.02. In contrast, 
the real PTF object distributions have extended tails beyond those 
ellipticity values. 

ellipticities in this region) . To limit our ellipticity search 
to stellar binaries, we measured the proper motion for 
each object in our sample from the USNO-Bl, 2MASS 
and (where available) SPSS measured position s (Monet 
et al.||2003[ [Skrutskie et~al][2006| [York 2000|). We re- 
quired at least a '2.ba proper motion detection for selec- 
tion for observation foUowup, with target brightness lim- 
its of 14 < mpi < 16.0, set by the PTF saturation limit 
and Robo-AO requirements. We selected 44 sources at 
a wide range of ellipticities from the resulting target list 
for follow-up observations with Robo-AO. 

4.1. Robo-AO 

Robo-AO is a visible and near-infrared laser guide star 
adaptive optics s ystem specifically eng ineered for 1-3 m 



class telescopes (Baranec et al. 2012). The Robo-AO 



system located at the Falomar bO-inch telescope com 
prises an ultraviolet Rayleigh laser guide star, an inte- 
grated adaptive optics and science camera system, and 
a robotic control system. The system currently incorpo- 
rates both an electron-multiplying CCD and an InGaAs 
infrared array camera for imaging. The robotic system 
is designed for high-efRciency observing, allowing large 
numbers of targets to be imaged rapidly. The diffraction- 
limited visible-light capability of Robo-AO allowed us to 
make a similar-wavelength comparison between the PTF- 
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Fig. 7.— Ellipticity distributions of the 100,000 PTF targets in 
the Robo-AO binary discovery target area. 

measured ellipticities and the observed binary contrast 
ratios. 

4.2. Binarity verification 

We obtained Robo-AO images of the 44 targets on the 
nights of July 16-18 and August 4-7 2012 (UT). The July 
targets, with RAs around 13 hours, were observed with 
60-second total exposure times in the i-band filter; the 
August targets (RAs of 22-23 hours) used a long-pass fil- 
ter with a 600 nm cut-on to obtain increased signal com- 
pared to a bandpass filter. We operated Robo-AO with- 
out tip-tilt correction, instead relying on post-facto shift- 
and-add processing o f the individual frames a nd used the 

|2012|) to perform 



pipeline described by Law et al 



(2009a 



the image alignment and co-addition. Ihe Robo-AO sys- 
tem provided diffraction-limited resolution (~ 0.1 arcsec 
at these visible wavelengths) for all observed targets, in 
1-2 arcsecond seeing conditions; the entire target list was 
observed in a total of only ^2.5 hours. 

The high-resolution images of our test sample (Figure 
[s]) confirm the expected variation of binarity with mea- 
sured ellipticity. Targets with small ellipticities (e<0.01) 
are all confirmed to be single stars. We found that five 
of the six targets with 0.01 < e < 0.02 are binaries; all 
the systems in that range are high-contrast and/or <1 
arcsecond separation. Targets with e > 0.02 are all con- 
firmed to be binaries, with separations increasing and 
contrast ratios decreasing as the ellipticities increase. A 
single target, PTF13.341, is detected as a binary with a 
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• 


PTF13.385 
e=0.0083 Vj 

* 


PTF23.583 ^ 
e=0.011 

# 


PTF22.443 
e=0.014 — ' 


PTF23.503 
e=0.014 U 

• 


PTF22.315 ^ 
• 


PTF13.341 A 
e=0.015 U 

« 


PTF13.341 A 
.e=0.017 \J 

• 


PTF22.255 
e=0.021 L/ 

« 


PTF23.593 
e=0.023 ^ 

' • 


PTF13.324 A 
e=0.025 VJ 


PTF22.302 
e=0.03 ^-J 


PTF22.463 ^ 
e=0.031 C-^ 


PTF23.512 (N 
e=0.032 Vj 

■ 


PTF14.591 
e=0.044 (-^ 

m 

* 


PTF23.553 (-^ 
-»* 


PTF22.303 , — . 
e=0.049 — ' 


PTF22.562 /-) 
e=0.052 L/ 


PTF13.315 A 
e=0.056 \J 

« 


PTF22.391 r\ 
e=0.058 \j 

s 


PTF22.551 ^ 
e=0.064 


PTF22.513 ^ 
e=0.069 

t 


PTF22.573 r\ 
e=0.072 \j 

• 

« 


PTF23.513 n\ 
e=0.074 V) 


PTF22.291 f\ 
e=0.091 V) 

• 

■ 


PTF23.570 
e=0.092 . U 

f 


PTF13.392 A 
e=0.11 \J 

t 


PTF22.243 
e=0.12 

» 


PTF22.591 ^ 
e=0.13 ~ 

• 


PTF13.225 . — , 
e=0.14 ^ 

• 


PTF22.594 A 
e=0.18 V) 
• 

• 


PTF13.294 pi 
e=0.22 ^ [J 


PTF13.22* A 
e=0.32 y 

• 



5 10 5 10 5 10 5 10 

X / arcsec 



Fig. 8. — Robo-AO adaptive-optics images of the binarity test targets. Each image is 10 by 10 arcseconds in size; East is towards the left 
and the images are rotated by 23.5° counter-clockwise with respect to a North-East axis. The images are shown with linear scaling with 
levels selected to best display the binarity (or lack thereof) of each of the targets. The ellipses indicate the predicted binary orientation 
based on the relative magnitudes of the ei and 62 parameters; the color indicates confidence level of the predicted companion (red: no 
companion; yellow: possible companion; green: very likely companion). The order of targets is the same as that of Table 1. 
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TABLE 1 

ROBO-AO TEST TARGETS, ORDERED BY INCREASING PTF-MEASURED ELLIPTICITY 





RA 


DEC 








separation 


r- 


band contrast 


Target Name 


(J2000) 


(J2000) 


ei 


62 


e 




(arcsec) 




(flux ratio) 


PTF13.262 


13:26:29.9 


+14:43:49 


-6.8e-05 


-4.3e-05 


8.1e-05 










Fib 14.544 


14:54:47.0 


+36:27:00 


n nnm a 

-0.00014 


n nnn "1 

0.00016 


n nnnoi 

0.00021 










Fib 13.2o4 


"1 O O A A t 

13:26:44.1 


1 T A OO rfcO 

+ 14:32:09 


n nnnoo 

0.00038 


n nnn o r' 

0.00025 


n n nn /I n 

0.00046 










Fib 22.385 


o o no 1 

22:38:56.1 


1 1 O A A A ~i 

+18:44:41 


-0.0005 


n nnn o "1 

0.00021 


n n nn K r' 

0.00055 










P it 23.593 


23:59:37.6 


1 o^ ^ o r" 

+27:48:55 


n nnnn"! 

0.00091 


0.00093 


0.0013 










Fib 22.021 


o o o "1 r o 

22:02:15.3 


+21:27:53 


-0.0015 


n nnn'^o 

0.00078 


0.0017 










Fib 23.443 


O O A A or "1 

23:44:35.1 


+28:00:39 


0.0019 


0.00092 


0.0021 










P 1 1 22.2S2 


OO OO OJ^ 

22:28:29.7 


1 "1 o n o 

+ 18:07:36 


-0.0021 


n nnn /i 

-0.00046 


n nnn o 

0.0022 










T~t rri T-\ o rzOA 

it zo. oo4 


OO cr o /in o 

23:58:40.8 


1 on n£? /I 

+29:06:46 


n nnno£? 

-0.00086 


0.002 


0.0022 










P iJ? 22.415 


OO A "I o o 

22:41:53.8 


+18:10:51 


-0.00076 


-0.0046 


0.0047 










P IJ? 22.343 


OO O A OO 

22:34:32.6 


1 1 o OO m 

+18:22:01 


n noo/^ 

0.0026 


n no j4 o 

0.0048 


0.0055 










nrriTTiOO cc(\ 

tr zz.ooU 


OO. c c .n 1 

22:&o:04.1 


+22:36:37 


0.0017 


-0.0053 


0.0056 










F lb 13.385 


1 *1 O O ^ o 

13:38:57.8 


1 A ^ 1 O O F 

+41:18:25 


n no ,( n 

-0.0049 


0.0067 


o noo o 

0.0083 










PTF23.583 


23:58:39.2 


+27:05:38 


0.01 


-0.0027 


0.011 





40 ± 0.04 





45 + 0.07 


PTF22.443 


22:44:30.6 


+18:34:25 


0.0074 


0.012 


0.014 





58 ± 0.04 





28 ± 0.03 


P 1 1 23.503 


23:50:34.0 


1 on n 'V on 

+29:07:39 


n nnnno 

0.00022 


-0.014 


0.014 


2 


01 ± 0.04 





n ■! 1 n f\ A 

04 ± 0.04 


P 1 1 22.315 


o o o "1 r* o o 

22:31:58.8 


1 "1 O A A A ^ 

+ 18:44:45 


0.012 


n nn '^n 

-0.0079 


0.014 


3 


78 + 0.04 





n '7 1 n n 

07 ± 0.06 


P 1 1 13.341 


"1 O O /I "1 O A 

13:34:18.4 


1 "1 o n o r* "1 

+ 13:08:51 


-0.012 


-0.01 


0.016 





20 + 0.04* 





on 1 n n 1 5k 

30 ± 0.01 




13:34:15.6 


1 TO in o o 

+ 13:10:38 


-0.015 


n nn '^o 

-0.0078 


0.017 












Zz.ZO.04. i 


1 1 n. c:'7. Qn 




-U.Ol / 


O.UZi 


1 


60 + 0.04 


u 


1 1 4- n no 


PTF23.593 


23:59:39.5 


+28:43:27 


-0.0018 


0.023 


0.023 





98 + 0.04 





17 + 0.02 


PTF13.324 


13:32:46.8 


+14:17:04 


-0.025 


-0.00032 


0.025 





96 + 0.04 





26 + 0.07 


PTF22.302 


22:30:25.1 


+18:22:13 


0.0041 


0.03 


0.030 





79 ± 0.04 





89 + 0.01 


PTF22.463 


22:46:36.0 


+18:59:59 


0.031 


-0.0037 


0.031 





73 ± 0.04 





52 + 0.01 


PTF23.512 


23:51:23.3 


+28:42:53 


-0.03 


0.011 


0.032 


2 


42 ± 0.04 





12 ± 0.02 


PTF14.591 


14:59:17.3 


+36:23:51 


0.039 


-0.02 


0.044 


3 


77 ± 0.04 





18 + 0.02 


PTF23.553 


23:55:35.6 


+29:01:03 


0.018 


0.044 


0.048 





84 ± 0.04 





43 + 0.02 


PTF22.303 


22:30:38.4 


+19:53:04 


0.032 


0.038 


0.049 





91 ± 0.04 





50 + 0.02 


PTF22.562 


22:56:24.9 


+23:15:57 


0.034 


-0.039 


0.052 





94 + 0.04 





91 + 0.01 


PTF13.315 


13:31:53.6 


+15:16:00 


-0.022 


-0.052 


0.056 


1 


64 + 0.04 





17 + 0.06 


PTF22.391 


22:39:14.8 


+19:17:56 


-0.047 


0.035 


0.058 


1 


07 + 0.04 





64 + 0.01 


PTF22.551 


22:55:16.4 


+22:30:49 


0.058 


-0.027 


0.064 


1 


01 ± 0.04 





55 + 0.01 


PTF22.513 


22:51:38.8 


+23:02:19 


0.024 


0.064 


0.069 


1 


92 + 0.04 





16 + 0.02 


PTF22.573 


22:57:37.7 


+22:58:29 


-0.069 


0.018 


0.072 


3 


67 + 0.04 





18 + 0.01 


PTF23.513 


23:51:30.5 


+28:26:14 


-0.06 


0.043 


0.074 


1 


25 + 0.04 





48 ± 0.01 


PTF22.291 


22:29:14.0 


+19:22:24 


-0.091 


0.0062 


0.091 


3 


81 ± 0.04 





14 + 0.03 


PTF23.570 


23:57:08.2 


+28:53:47 


0.018 


-0.09 


0.092 


2 


75 ± 0.04 





12 + 0.03 


PTF13.392 


13:39:21.8 


+41:25:23 


-0.078 


-0.075 


0.11 


1 


61 ± 0.04 





57 + 0.01 


PTF22.243 


22:24:37.1 


+19:55:23 


0.046 


0.11 


0.12 


4 


49 ± 0.04 





12 ± 0.05 


PTF22.591 


22:59:16.7 


+23:18:27 


0.12 


-0.061 


0.13 


2 


74 ± 0.04 





17 ± 0.03 


PTF13.225 


13:22:59.8 


+41:11:44 


0.12 


0.072 


0.14 


3 


04 ± 0.04 





17 + 0.03 


PTF22.594 


22:59:43.1 


+23:53:08 


-0.18 


0.012 


0.18 


2 


05 + 0.04 





62 + 0.01 


PTF13.294 


13:29:49.7 


+41:06:13 


-0.028 


-0.22 


0.22 


2 


29 + 0.04 





33 + 0.07 


PTF13.224 


13:22:41.6 


+14:00:14 


-0.1 


-0.3 


0.32 


3 


82 + 0.04 





43 + 0.02 



* — companion not associated with PTF ellipticity 
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faint, close companion by Robo-AO, but with a different 
position angle and closer separation than that measured 
from PTF data; this is likely to be a close companion not 
measured by our PTF data and we do not include it in 
further analysis. 

The crossover point between single stars and high- 
confidence binaries occurs at the ellipticities predicted 
by our artificial point source simulations (Section [3|, al- 
though the >80% binarity fraction of the targets with 
ellipticities between 0.01 and 0.02, compared to the zero 
binarity fraction at lower ellipticities, suggests that a 
less conservative limit could be set for many science pro- 
grams. 

5. CHARACTERIZING DETECTED BINARIES: POSITION 
ANGLES, SEPARATIONS AND CONTRAST RATIOS 

Beyond simple detection of binarity the measured el- 
lipticity vector allows us to directly measure the binary's 
position angle. In a multiple-epoch survey, the variation 
in ellipticity with seeing enable us to constrain both the 
binary's separation and its contrast ratio. 

5.1. Position angles 

The detected Robo-AO binaries show that our ellip- 
ticity measurements can predict the orientation of the 
binary from the relative magnitudes of the ei and 62 
parameters (see the blue ellipses in Figure [S]) . The ei 
parameter varies as cos(20) and the 62 parameter varies 
as sin(20), where 9 is the angle that the vector connect- 
ing the binaries makes with the horizontal. This makes 
it possible to predict 9 via 



tan 



(13) 



In our test dataset, the median difference between the 
PTF-predicted orientation and the observed Robo-AO 
orientation is 2.1°. 

5.2. Separations and contrast ratios 

The ellipticity of a close binary increases when the flux 
ratio of the stars in the binary system increases, or al- 
ternatively when the separation between the binaries in- 
creases (Figure [5| . Because of this degeneracy between 
separation and contrast, it is not possible to estimate the 
binary parameters by examining solely the magnitude of 
the ellipticity. However, we can differentiate between the 
two cases based on how the ellipticity varies with seeing. 

When the members of a binary are very close together, 
the PSF of their combined light will only be noticeably 
elliptical for smaller FWHM values. As the FWHM in- 
creases, the PSF will quickly start to resemble a circular 
PSF. On the other hand, when the two members of a bi- 
nary are separated by a large angular distance, their com- 
bined PSF will remain elliptical for much larger FWHMs 
than the close binaries. Thus, if the magnitude of the el- 
lipticity of a binary is relatively constant with seeing, 
the binary has a wide separation, and if the ellipticity 
decreases quickly with seeing, the binary members must 
be close together. 

For example, two of the binaries from Figure [8] 
PTF22.573 and PTF23.513, have nearly identical ellip- 
ticities at the chosen reference FWHM of 2.5 arcseconds 
(as measured from the PTF images). The high-resolution 



m 10- 
>- 

15 



PTF22.573 



PTF23.513 




2.4 2.6 
seeing (arcseconds) 

Fig. 9. — Comparison of the seeing dependence of the eUipticity 
of two binary systems which have a similar measured ellipticity at 
a common reference FWHM. The images of the two binary sys- 
tems are from the PTF survey. The ellipticity of the binary whose 
members are closer together is much more sensitive to the FWHM, 
while the ellipticity of the binary whose members are further apart 
is almost constant with changes in FWHM. This difference in see- 
ing dependence allows us to resolve the seeing/contrast degeneracy 
of the ellipticity. 

Robo-AO images reveal that these two binaries differ in 
contrast and separation - PTF22.573 has greater angu- 
lar separation and lower contrast ratio than PTF23.513. 
Figure |9] shows the seeing dependence of the ellipticity 
of both of these binaries. It is apparent than the closer, 
higher contrast binary has an ellipticity which is much 
more sensitive to the FWHM, as expected. This differ- 
ence in response to changes in seeing allows us to resolve 
the degeneracy in the ellipticity measurement between 
separation and contrast. 

For the binaries we found with Robo-AO, we examined 
how the measured ellip ticity in the PTF images varied 
with seeing (Figure 



10 1. We derived a simple linear fit of 



the ellipticity as a function of seeing for each of the bi- 
naries. We evaluated these fits at a reference FWHM of 
2.5 arcseconds, and we parametrized the response of the 
ellipticity to seeing with the fractional decrease in seeing 
per arcsecond (Ae / e, where Ae is the magnitude of the 
slope of the linear fit of ellipticity with respect to seeing, 
and e is the ellipticity measured at the reference radius 
of 2.5 arcseconds). As can be seen in Figure 10 the bi- 



naries with widest angular separations have close to zero 
variance of ellipticity with seeing. Using the measured 
relation to predict the Robo-AO measured separation on 
the basis of PTF data alone confirms BinaryFinder can 
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1.0 2.0 3.0 4.0 
separation (arcsec) 



5.0 




1.0 2.0 3.0 4.0 
PTF-estimated separation (arcsec) 



00 0.02 0.04 0.06 0.08 0.10 0.12 
PTF-e / PTF estimated separation 



Fig. 10. — Left: the change in eUipticity with seeing (as a fraction of the magnitude of ellipticity) for the binaries we imaged with Robo- 
AO. For larger separations, the ellipticity becomes less responsive to the FWHM. Middle: comparison between the separations estimated 
from PTF data alone and the Robo-AO measured separations; the correlation shows separation measurement at ~25% precision. Right: 
comparison between the PTF ellipticity divided by PTF separation, and the Robo-AO measured flux ratio. Although the correlation 
is weak, there is a significant excess of points towards the lower left, suggesting that the contrast ratio of the binaries can be weakly 
constrained using PTF data alone. 



measure the binary separations to ^-^25% precision. 

The measured ehipticity is a function of flux ratio and 
separation only, and so with a measured separation the 
flux ratio of the binary can also be estimated. Also in 
Figure [TOj we use the estimated separations to constrain 
the flux ratios of the systems, assuming a simple linear 
relationship between separation and ellipticity indepen- 
dent of contrast ratio (which is approximately correct; 
see Sectionjsl) . The resulting distribution is (weakly) cor- 
related withthe Robo-AO measured flux ratios; although 
there are outlier points, a low value of ellipticity-over- 
separation strongly suggests a high flux ratio between 
the components of the binary systems. 

6. DISCUSSION, APPLICATIONS AND CONCLUSIONS 

6.1. False positives 

False positives can be generated in our algorithm by in- 
strumental effects that mimic PSF ellipticity, or extended 
astrophysical objects which are not physically- associated 
binaries. 

The careful attention to input data quality and PSF 
anisotropy removal results in a low level of instrumental 
false positives. Of the 44 targets observed with Robo- 
AO, 25 had ellipticities greater than our blend detec- 
tion criterion of 0.02 ellipticity. All 25 of these targets 
were confirmed to be binary systems, suggesting a false- 
positive rate of less than 5% level. 

Proper motion sample cuts can ensure that the Bina- 
ryFinder targets are relatively nearby stars rather than 
galaxies and other extended extragalactic objects. We 
note, however, that if the proper motion cuts are relaxed, 
the technique detailed in this paper can detect and mea- 
sure the characteristics of any objects that lead to an 
elliptically-shaped image. Beyond binary detection, the 
algorithm could thus be used to search for and statisti- 
cally characterize faint barely-extended galaxies, strong 
gravitational lenses and other extragalactic objects. 

6.2. Very-large-sample binary detection 

We have demonstrated that BinaryFinder can detect 
binaries down to «l/5 of the seeing limit in wide-field 
survey data, and directly measure their position angles, 
separations and contrast ratios. When applied to the 
synoptic sky survey data currently being collected in 
great quantities by synoptic sky surveys like PTF and 



Pan-Starrs, we will be able to perform a search for bi- 
nary systems over very large sky areas. 

At the time of writing PTF has imaged 1200 fields with 
sufhcient numbers of epochs for full BinaryFinder binary 
detection and measurement of position angles and sepa- 
rations. Averaging across the sky, each PTF field con- 
tains «10,0G0 objects which are bright enough (m^ < 18) 
for the detection of companions at separations down to 
~l/5 of the seeing limit. In typical fields, we find that 
approximately 35% of those objects have high enough 
proper motions to confirm that they are nearby stars, 
and wll% of those are detected as < 2 arcsecond separa- 
tion binaries by our algorithm (either binaries or blended 
unassociated foreground/background stars). Estimating 
the background blend probability from the stellar den- 
sity in our images, we estimate that only «1% of those 
targets will be background blends. 

In total, a full PTF-datasct BinaryFinder search could 
target ^12,000,000 stars; we would expect to detect 
^450,000 binary systems with detectable proper motions 
and a high confidence of physical association. A binary 
sample of that size will allow unprecedentedly detailed 
statistical analysis of the binarity fraction as a function 
of stellar mass, age, metallicity and galactic population. 

PTF and similar surveys operate in multiple filters. 
Multi-color imaging ellipticity measurements, if avail- 
able, would scale according to the relative colors of the bi- 
nary components, and could be used to reduce the false- 
positive probability or further constrain the nature of 
detected systems. 

6.3. Orbital Measurements for an Extremely Large 
Binary Sample 

Many current and planned synoptic sky surveys image 
the same areas of sky repeatedly over many years. With 
the ability to directly measure position angle, separation 
and contrast ratio of individual binaries, in many cases 
sufficient data will be collected to follow the motion of 
binary systems with time. Since this can be done over 
large areas and very large samples, BinaryFinder could 
be used to obtain orbital measurements (and thus mass 
constraints) for a very large sample of binaries. 

6.4. Summary and conclusions 
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Wc have shown that BinaryFindcr is capable of de- 
tecting and characterizing close binary systems in wide- 
field synoptic survey data. Using Robo-AO, we have con- 
firmed the PTF direct detection and separation measure- 
ment of binaries as small as 0.4 arcsec, or less than 1/5 of 
the PTF FWHM. Images from the PTF camera have a 
pixel sampling of 1.01 arcseconds and a median FWHM 
of «2.0 arcseconds. Because our algorithm requires well- 
sampled images, we restricted its operation to PTF im- 
ages with 2-3 arcsec FWHMs. New synoptic sky surveys 
with smaller pixels and smaller PSFs will enable the de- 
tection of much closer binary systems; for well-sampled 
images the detection limit scales with the PSF FWHM 
(assuming similar systematic noise performance). Using 
BinaryFinder, a survey with improved sampling such as 
LSST is likely to be capable of the detection, character- 
ization and orbital motion measurement of millions of 
multiple systems with separations as close as 0.1 arcsec- 
onds. 
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